Probing light WIMPs with directional detection experiments 
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The CoGeNT and CRESST WIMP direct detection experiments have recently observed excesses 
of nuclear recoil events, while the DAM A/LIBRA experiment has a long standing annual modulation 
signal. It has been suggested that these excesses may be due to light mass, m x ~ 5—10 GeV, WIMPs. 
The Earth's motion with respect to the Galactic rest frame leads to a directional dependence in 
the WIMP scattering rate, providing a powerful signal of the Galactic origin of any recoil excess. 
We investigate whether direct detection experiments with directional sensitivity have the potential 
to observe this anisotropic scattering rate with the elastically scattering light WIMPs proposed to 
explain the observed excesses. We find that the number of recoils required to detect an anisotropic 
signal from light WIMPs at 5a significance varies from 7 to more than 190 over the set of target nuclei 
and energy thresholds expected for directional detectors. Smaller numbers arise from configurations 
where the detector is only sensitive to recoils from the highest speed, and hence most anisotropic, 
WIMPs. However, the event rate above threshold is very small in these cases, leading to the need 
for large experimental exposures to accumulate even a small number of events. To account for this 
sensitivity to the tail of the WIMP velocity distribution, whose shape is not well known, we consider 
two exemplar halo models spanning the range of possibilities. We also note that for an accurate 
calculation the Earth's orbital speed must be averaged over. We find that the exposures required 
to detect 10 GeV WIMPs at a WIMP- proton cross-section of 10 -4 pb are of order 10 3 kgday for a 
20 keV energy threshold, within reach of planned directional detectors. Lower WIMP masses require 
higher exposures and/or lower energy thresholds for detection. 

PACS numbers: 95.35.+d 



I. INTRODUCTION 

Direct detection experiments aim to detect dark mat- 
ter in the form of Weakly Interacting Massive Parti- 
cles (WIMPs) via the nuclear recoils which occur when 
WIMPs scatter off target nuclei pj. The sensitivity of 
these experiments has increased rapidly over the last few 
years, and they are probing the regions of WIMP mass- 
cross-section parameter space populated by the lightest 
neutralino in Supersymmetric extensions of the standard 
model (see e.g. Ref. [2[). 

Event rate excesses and annual modulations in various 
direct detection experiments have prompted recent inter- 
est in light WIMPs. The DAMA (now DAM A/LIBRA) 
collaboration have, for more than a decade, observed an 
annual modulation of the event rate in their Nal crys- 
tals Q. This annual modulation is consistent with light 
(m x ~ 5 - 10 GeV) WIMPs scattering off Na d i]. 
The CoGeNT experiment, after allowing for backgrounds 
with an exponential plus constant energy spectrum, find 
an excess of low energy events which is consistent with 
WIMPs with mass m x » 7 - 11 GeV 6]. With a larger 
data set they have observed a 2.8 a annual modula- 
tion Q, with period and phase broadly consistent with 
the expectation for WIMPs Q. The CRESST experi- 
ment has observed an excess of events in their CaWC>4 
crystals above expectations from backgrounds [9]. The 
excess is compatible with either WIMPs of mass m x ~ 25 
GeV scattering off tungsten predominantly, or lighter, 
m x ~ 10 GeV, WIMPs scattering off oxygen and calcium. 
It appears that it is not possible to explain all of these 



signals in terms of a single conventional elastic-scattering 
WIMP, esp ecially when the exclusion limits from the 
CDMS Hflf, XENON10 [uj and XENON100 [II experi- 
ments and the CRESST commissioning da ta J13II data are 
taken into account 0-GJ] (see also Refs. dMl). None 
the less it is still possible that some subset of the putative 
signals are due to elastic scattering light WIMPs. 

The deployment of a Nal detector at the South Pole 
has been proposed to directly test the DAMA annual 
modulation signal [l{|. The direction dependence of the 
scattering rate [2(| provides another potentially powerful 
way of testing whether the observed excesses and annual 
modulations are due to elastic scattering light WIMPs. 
The amplitude of the directional signal is far larger than 
that of the annual modulation and hence the anisotropy 
of the WIMP induced nuclear recoils could be confirmed 
with a relatively small number of events [2lL [22} . Fur- 
thermore the angular dependence of the recoils (in par- 
ticular the peak recoil rate in the direction opposite to 
the motion of the solar system, or for low energy recoils 
a ring around this direction [231 ]) is extremely unlikely 
to be mimicked by backgrounds, and would allow unam- 
biguous detection of WIMPs 0, [25| . In this paper we 
investigate whether current and near future directional 
detectors would be able to detect elastic scattering light 
WIMPs. 



II. MODELLING 

We use the same statistical techniques and methods for 
calculating the directional nuclear recoil spectrum as in 
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Refs. [22], |26|, |27j . We briefly summarise these procedures 
here, for further details see these references. 



A. Detector 

Most of the directional detectors currently under de- 
velopment (see Refs. [2^,[2^] for reviews) are low pressure 
gas time projection chambers (TPCs), e.g. DMTPC (H^, 
DRIFT [U, MIMAC [H and NEWAGE Various 
gases have been considered, including CF4, CS2 and 3 Hc. 
We therefore consider all four of these target nuclei: 3 He, 
C, FandS. 

Detailed calculations of the nuclear recoil track recon- 
struction are not available for all of these targets (see 
Ref. [lil for a detailed study of the reconstruction of 
simulated tracks for a MIMAC-like detector). Therefore 
we assume that the recoil directions are reconstructed 
perfectly in 3d. This is an optimistic assumption, there- 
fore our results provide a lower limit on the number of 
events and exposure required by a real TPC based detec- 
tor. Finite angular resolution does not significantly affect 
the number of events required to detect the anisotropic 
WIMP signal, provided it is not worse than of order tens 
of degrees (H,!!!, H3|- 2d read-out would, however, sig- 
nificantly degrade the detector capability [22|, [26|, [27|, [35| . 
We consider both vectorial and axial data i.e. where the 
senses of the recoils +x and — x are either measured for 
all recoil events or no events . Sense discrimination is 
a major challenge for directional detectors. As discussed 
in detail in Ref. [Hj], while the shape and charge distri- 
bution of nuclear recoil tracks are expected to be asym- 
metric, measuring these asymmetries with high efficiency 
is difficult in practice. 

We consider four benchmark energy thresholds: E t h — 
5, 10, 15 and 20 keV for each target. Note that these 
are directional energy thresholds. It is harder to mea- 
sure the direction of a recoil than to simply detect it 
therefore, for a given experiment, the directional energy 
threshold is usually larger than the threshold for simply 
detecting recoils. The high energy recoils are the most 
anisotropic [2(|, therefore for heavy WIMPs a low en- 
ergy threshold is not essential for directional detection. 
For instance for WIMPs with m x = 100 GeV and a S 
target, reducing the energy threshold below 20 keV does 
not significantly reduce the exposure required to reject 
isotropy |27j. However for light WIMPs the differential 
event rate decreases rapidly with increasing energy, and 
a low threshold is crucial to obtain a non- negligible event 
rate. We discuss the viability of sufficiently low energy 
thresholds in Sec. IIIII 



B. WIMP masses and cross-sections 



We consider three benchmark light WIMP masses, 
m x = 5,7.5 and 10 GeV spanning the range of masses 
where the observed nuclear recoil excesses might be con- 
sistent with exclusion limits from other experiments. We 
fix the elastic scattering cross-section on the proton to 
(T p = 10~ 4 pb. It is straight-forward to scale our re- 
sults to other values for the cross-section. The number of 
events required to detect anisotropy, N- lso , is independent 
of the cross-section, while the corresponding exposure, £, 
is calculated as 
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a pP R(> Eth) ' 



(1) 



1 For studies of the effects of statistical sense determination see 
Refs. [H, [Hi. 



where R(> E t h) is the total event rate (i.e. the integral 
of the differential event rate) above the energy threshold 
normalised to unit cross-section and local WIMP density. 
Therefore the exposures required can be simply scaled for 
other cross-sections and local WIMP densities. 



C. WIMP velocity distribution 

The detailed angular dependence of the recoil rate de- 
pends on the exact form of the WIMP velocity distribu- 
tion [H [H, However, if the WIMP velocity distri- 
bution is dominated by a smooth component the main 
feature of the recoil distribution (namely the rear-front 
asymmetry) is robust (see e.g. Ref. [22])). The number 
of events required to detect anisotropy depends relativity 
weakly on the WIMP speed distribution [22j ]. However 
the event rate above the energy threshold, and hence the 
exposure required to detect anisotropy, depends more sig- 
nificantly on the WIMP speed distribution. 

Usually the dominant uncertainty in the (time and di- 
rection averaged) differential event rate comes from the 
~ 10% uncertainty [4(| in the value of the local cir- 
cular speed, and hence the WIMP velocity dispersion 
(e.g. Refs. [SWI^ y The velocity dispersion and WIMP 
mass have somewhat degenerate effects on the differential 
event rate. For instance if the velocity dispersion is in- 
creased, then there are more WIMPs with higher speeds, 
however the energy spectra of the resulting nuclear re- 
coils can remain the same if the WIMP mass is decreased. 
Therefore the range of WIMP masses corresponding to 
a particular observed energy spectrum excess moves to 
lower masses if the velocity dispersion is increased (see 
e.g. Ref. 17] for the specific case of CoGeNT). Vary- 
ing the velocity dispersion has a qualitatively similar ef- 
fect on the values of the WIMP mass consistent with 
the DAMA annual modulation [45[ . Consequently, while 
varying the circular speed affects the values of the WIMP 
mass corresponding to, or excluded by, the various data 
sets, it does not significantly affect their compatibility. 
Therefore we fix the local circular speed to its standard 
value, v c — 220kms~ 1 , consistent with our benchmark 
WIMP masses m x = 5, 7.5 and 10 GeV. 
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FIG. 1: The minimum WIMP speed, v m i n , which can cause 
a recoil of energy E for 3 He (solid), C (dotted), F (short 
dashed) and S (long dashed) for (top and bottom curves in 
each case respectively) m x = 5 and 10 GeV. The horizontal 
dot-dashed lines show the maximum WIMP speeds in the lab, 
^max w v ^ _|_ 248km s _1 , corresponding to the 90% upper 
and lower confidence limits on the escape speed from RAVE, 
v™ x = 608 km s" 1 and v™ n = 498 km s" 1 . 



Since the differential event rate involves an average 
over the WIMP speed distribution it is usually relatively 
weakly sensitive to the detailed shape of the speed dis- 
tribution (e.g. Refs. [12, Hi|). This is not necessarily 
the case, however, for experiments that are only sen- 
sitive to the high speed tail of the distribution (e.g. 
Refs. [II III Ell). The minimum WIMP speed which 
can cause a recoil of energy E, v m - m , is given by 
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where m A is the mass of the target nuclei and P x a is the 
reduced mass of the WIMP-target nucleus system. 

Particles with speed greater than the local escape 
speed, v csc = ^/2|(f>(_R ) where $ is the potential and Rq 
the Solar radius, will not be gravitationally bound to the 
Milky Way 2 . The RAVE survey found that the escape 
speed lies in the range 498 km s" 1 < v csc < 608 km s -1 
at 90% confidence, with a median likelihood of v csc — 
544 km s" 1 [H. The 

maximum WIMP speed in the lab 
frame is v csc + v\ a h(t), where i>i a t>(i) is the speed of the 
Earth with respect to the Galactic rest frame. This is 
made up of three components: the motion of the Local 
Standard of Rest (LSR), vlsr = (0, u c ,0) in Galactic 
coordinates, the Sun's peculiar motion with respect to 
the LSR, v|(ll.l, 12.2, 7.3) kms" 1 [H, and the Earth's 



2 Dark matter halos do contain unbound particles, however the 
fraction of such particles at the Solar radius is small O . 



orbit about the Sun, v° rb (i). It has a maximum value 
at t ks 153 days (on June 2nd) of v"^ x = Wiab(*o) ~ 
248 kms" 1 . 

Fig. [T] shows w m in as a function of E for m x = 5 and 
10 GeV for 3 He, C, S and F target nuclei. The hor- 
izontal lines show the maximum WIMP speeds in the 
lab, «™ ax = ti osc + f™b X i corresponding to the 90% up- 
per and lower confidence limits on the escape speed from 



RAVE, 



608 kms and u' 1 



498 kms" 



For 



light WIMPs, unless the target nuclei are light and the 
threshold energy low, the minimum speed corresponding 
to the threshold energy lies in the tail of the speed dis- 
tribution, and in some cases beyond the cut-off due to 
the Galactic escape speed. The expected event rates are 
therefore very sensitive to the value of the escape speed 
and the shape of the high speed tail of the distribution. 
The standard halo model, an isotropic, isothermal 



sphere with density profile p(r) oc 



is formally in- 



finite. Hence its Maxwellian velocity distribution, 



./<vi= V,. M ,(-^£ 
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where a — y'3/2v c and N is a normalisation constant, 
extends to infinity too. This is usually addressed by trun- 
cating the velocity distribution by hand, either sharply 
or exponentially, at the escape speed. 

Numerical simulations find velocity distributions 
with less high sp eed particles than the standard 
Maxwellian [51-53]. Lisanti et al. [47j have presented 
an ansatz for the velocity distribution which reproduces 
this behaviour: 
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eK S c-|v|). (4) 



The parameter k is related to the outer slope of the den- 
sity profile, 7, (p(r) oc r~ 7 for large r), by k = 7 — 3/2 
for 7 > 3 [Hl j. 

We consider 2 forms for /(v) chosen to have high speed 
tails which roughly span the plausible range. Firstly, 
a Maxwellian distribution, eq. with a sharp cut-off 
(/(v) = for |v| > v™ x ) at the upper limit on the es- 
cape speed from RAVE, v™* = 608 kms -1 , which pro- 
vides a large tail event rate. Secondly, a Lisanti et al. 
f(v), eq. flU, with v csc = w™ n = 498kms _1 to pro- 
vide a small tail event rate. We fix k — 1.5, corre- 
sponding to an outer density profile slope 7 = 3 and 
vq = v c — 220 kms . In both cases we fix the local 
density to be p = 0.3 GeV cm -3 . Scaling the event rates 
and exposures to other densities is straight-forward. 

Fig. [2] shows the differential event rates for each of 
the WIMP masses and targets we consider calculated us- 
ing the Maxwellian distribution with a sharp cut-off at 



= 608 kms" 



The much smaller event rates for 



He are largely due to the A 2 factor in the event rate for 
spin-independent scattering. For heavier targets the dif- 
ferential event rate in the E — > limit is substantially 
larger, however it decreases rapidly with increasing E, 
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FIG. 2: The differential event rates, assuming a Maxwellian 
speed distribution with a sharp cut-off at v™J x = 608 kms -1 , 
a cross-section <r p = 10 -4 pb and a local WIMP density p = 
0.3GeVcm" 3 , for He (solid), C (dotted), F (short dashed) 
and S (long dashed) for (from top to bottom at E — OkeV) 
m x = 5, 7.5 and 10 GeV. 

and the lighter the WIMP the more rapid the decrease. 
Furthermore for relatively small energies u m in exceeds the 
maximum WIMP speed in the lab frame and hence the 
differential event rate is zero. 

Fig. [3] shows the ratio of the speed integral, 




for Lisanti et al.'s f(v) in eq. @ with w™ n = 498 kms -1 
to that for the Maxwellian distribution with a sharp cut- 
off at v"lc* — 608 kms -1 . It also shows the velocity in- 
tegral ratios for these two models if the Earth's orbit is 
neglected. For v m i n < 0(300 kms -1 ) the difference in 
the speed integrals for the two speed distributions is rel- 
atively small, less than 10%. As w m i n is increased, so 
that only the tail of the speed distribution is included 
in the speed integral, the difference becomes substan- 
tially larger, reaching roughly an order of magnitude for 
Vmin ~ 650 kms -1 . Neglecting the Earth's orbit has a 
significant (> 10%) effect, for i> m i n within ~ 100 kms -1 
of the maximum WIMP speed in the lab frame w™ ax . 
Therefore the Earth's orbital speed must be included, 
and averaged over, for an accurate calculation of the 
event rate. 

Simulated dark matter halos have velocity distribu- 
tions which contain features at high speeds HH . More 
specifically there are fairly broad features, which are sim- 
ilar at different positions within a single halo, but which 
vary from halo to halo, and are hence thought to be a relic 
of the formation history of the halo [H, HH . Ref . [53[ 
also finds narrow features in some locations, correspond- 
ing to tidal streams, while Ref. 56] finds that the dark 
matter streams from the Sagittarius dwarf are signifi- 
cantly more extended than the stellar streams, and the 



FIG. 3: The speed integral, g(v m i n ), relative to that for 
the Maxwellian distribution with a sharp cut-off at f™c X = 
608 kms -1 , including averaging over the Earth's orbit. The 
solid line is for Lisanti et al.'s /(«), eq. Q, with w™c" = 
498 kms -1 , including the Earth's orbit. The dotted and 
dashed lines are for the Maxwellian and Lisanti et al. dis- 
tributions respectively, neglecting the Earth's orbit. 



leading dark matter stream may pass through the Solar 
neighbourhood. The detailed shape of the high speed tail 
of f(v) would affect the interpretation of the CoGeNT, 
CRESST and DAMA data, in particular the values of 
the WIMP mass extracted. Since no detailed study of 
this has been carried out to date 3 , we do not include the 
high speed features in our analysis. We defer a general 
investigation of the directional event rate produced by 
simulation velocity distributions to future work (58j . 



D. Statistical tests 

We follow the statistical procedures described in detail 
in Refs. [22|, |27j. We use the Rayleigh- Watson statistic, 
which uses the mean resultant length of the recoil direc- 
tion vectors. We also use the Bingham statistic which, 
unlike the Rayleigh- Watson statistic, can be used with 
axial data (where the senses of the nuclear recoils are 
not measured). These statistics are described in more 
detail in Appendices [A] and |B] 

For each WIMP mass, target nuclei and energy thresh- 
old combination we calculate the probability distribution 
of each statistic for WIMP induced recoils and also for 
the null hypothesis of isotropic recoils. We use these 
distributions to calculate the rejection and acceptance 



3 Ref. [l5l | includes tidal streams with specific properties chosen to 
reproduce the energy dependence of the amplitude of the annual 
modulation measured by CoGeNT, while Ref. [57j studies the 
effects of the Sagittarius tidal stream. 
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factors, R and A. The rejection factor gives the con- 
fidence level with which the null hypothesis can be re- 
jected given a particular value of the test statistic, while 
the acceptance factor is the probability of measuring a 
larger value of the test statistic if the alternative hypoth- 
esis is true. We then find the number of events required 
for A = R = 0.95,0.99730 and 0.999 999 427 i.e. to re- 
ject isotropy at this confidence level in this percentage 
of experiments (see Refs. [H, H3 for further discussion). 
The later two confidence levels correspond, for a gaussian 
distribution, to three and five sigma respectively. Since 
our aim is to examine whether directional detection ex- 
periments could detect light WIMPs we will focus on the 
case A = R = 0.999 999 427, corresponding to the 'five 
sigma' result conventionally required for discovery. 
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III. RESULTS AND DISCUSSION 

For each of the combinations of WIMP mass, target 
nuclei, energy threshold and confidence level discussed 
in Sec. [TT] we calculate the number of events required 
to reject isotropy, -/Vj so , for vector and axial data (us- 
ing the Rayleigh- Watson and Bingham statistics respec- 
tively). The number of events required for a 5<r detection 
with the Rayleigh- Watson statistic varies from 7 to 58. 
For fixed WIMP and target mass, iV; so decreases with in- 
creasing energy threshold, E t h, (since the recoils caused 
by high speed WIMPs in the tail of the distribution are 
more anisotropic), until the minimum speed required to 
cause a recoil of energy Eth, v m - m (Eth), exceeds the max- 
imum WIMP speed in the lab frame, i>™ ax . At this point 
the event rate is zero and no events can be detected. 
The same trend occurs for decreasing WIMP mass (with 
threshold energy and target mass fixed). For fixed thresh- 
old energy, iVj so decreases with increasing target mass for 
m x = 5 GeV, however for m x — 7.5 and 10 GeV, Ai so 
increases as the target mass number is increased from 
A = 3 to 12 before decreasing as A is increased further. 
This is due to the variation of v m i n with target and WIMP 
mass shown in Fig. [TJ 

For the Bingham statistic, which can be used with ax- 
ial data, the number of events required for a 5<r detec- 
tion varies from 9 to more than 190 4 . Both the num- 
ber of events and its increase, relative to the number re- 
quired for the Rayleigh- Watson statistic, is smallest for 
the configurations which are only sensitive to the highly 
anisotropic recoils from high speed WIMPs in the tail of 
the speed distribution. 

If an experiment is only sensitive to high speed 
WIMPs, fewer events are required to reject isotropy, how- 
ever, the reduced event rate means that the exposure re- 



4 In a small number of cases, where a large number of events are 
required, we have only been able to place a lower limit on N- lao 
due to computational time limitations. 



FIG. 4: The exposure required for a 5cr rejection of isotropy as 
a function of energy threshold, -Eth, for 3 He (top left), C (top 
right), F (bottom left) and S (bottom right) for m x = 5,7.5 
and 10 GeV (crosses, stars and triangles respectively) for the 
Maxwellian f(v) with a sharp cut-off at = 608 km s . 

Where a symbol is not displayed, the minimum WIMP speed 
corresponding to the energy threshold, v m i n (Eth), for this 
WIMP and target mass combination exceeds the maximum 
WIMP speed in the lab frame, w™ ax and hence the event rate 
is zero. 

quired to accumulate these events will be larger. We 
therefore use eq. ([T]) to calculate the exposure, £, (in 
kg day) required to accumulate the required number of 
events for each case. As illustrated in Fig. [2 for light 
WIMPs the differential event rate decreases rapidly with 
increasing energy, and therefore the event rate above the 
energy threshold, R(> E t h), plays a crucial role in deter- 
mining the exposure required. 

The exposure required to reject isotropy at 5er us- 
ing the Rayleigh- Watson statistic assuming a Maxwellian 
f(v) with a sharp cut-off at f™ ax = 608 km s -1 is shown 
for each configuration in Fig. HI While jVj so varies by 
less than an order of magnitude, because of the large 
variation in R(> -Eth), the exposure varies by more than 
five orders of magnitude. Due to the rapid decrease of 
R(> E t h), the exposure increases sharply with increas- 
ing E t h for each WIMP and target nuclei mass combi- 
nation. The factor by which the exposure increases, in- 
creases with both decreasing WIMP mass and increasing 
target nuclei mass (i.e. as the minimum WIMP speed to 
which the experiment is sensitive is increased). Eventu- 
ally the minimum WIMP speed corresponding to the en- 
ergy threshold, v m - m (Eth), exceeds the maximum WIMP 
speed in the lab frame, v™ ax , and the event rate is zero 
and WIMPs of this mass can not be detected. 

Of the halo models considered, the Maxwellian f(v) 
with a sharp cut-off at v™ ax = 608 km s^ 1 has the largest 
tail event rate, and hence the smallest exposures. In 
Fig. [5]we show the exposures for a S target for the Lisanti 
et al. f(v), eq. Q, with w™ n = 498 km s" 1 as well. 
When v m i n (E t h) is much smaller than ?j™ ax the expo- 
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FIG. 5: As Fig. Ubut for a S target only comparing the expo- 
sures required for the Lisanti et al. f(v), eq. ©, with v™ 1 ™ = 
498 km s -1 (upper symbols) with those for the Maxwellian 
f(v) with a sharp cut-off at i>2c X = 608 km s -1 (lower sym- 
bols). For m x = 7.5 keV and E t h = 5keV the exposures 
required for the two speed distributions are the same, and 
hence only one symbol is visible. 



FIG. 6: As Fig.[S]but comparing the exposures required using 
the Bingham statistic (upper symbols) with those required 
using the Rayleigh statistic (lower symbols), in both cases for 
a S target using a Maxwellian f(v) with a sharp cut-off at 
= 608 km s" 1 . For m x = lOGeV and E th = 5keV we 
have only been able to place a lower bound on the number 
of events, and hence exposure, required with the Bingham 
statistic. 



sure required is fairly modest, S ~ 10 — 100 kg day and 
the event rates, and hence exposures, for the two speed 
distributions are very similar. However as v m i n (E t h) ap- 
proaches v™ ax the exposures required become large and 
the differences between the two speed distributions be- 
come significant. For instance for E t h = 20keV, WIMPs 
with m x = 10 GeV and a Maxwellian distribution with 
v clc X = 608 km s -1 anisotropy could be detected with an 
exposure of 1700 kg day, however if the WIMPs have the 
Lisanti et al. f(v) with w™ n = 498 km s" 1 the event rate 
is zero and they can not be detected. The trends for the 
other target nuclei are similar. 

In Fig. [5] we compare the exposures required to reject 
isotropy with axial data using the Bingham statistic with 
those for vectorial data using the Rayleigh- Watson statis- 
tic, for a S target and a Maxwellian f(v) with a sharp 
cut-off at vZc* = 608 kms -1 . Since R(> E th ) for each 
configuration doesn't change, the variations in the expo- 
sure are driven entirely by the variations in 7V; so discussed 
above. Therefore the increase in the exposure, relative to 
that required for the Rayleigh- Watson statistic, is small- 
est for the cases which are only sensitive to the highly 
anisotropic recoils from high speed WIMPs in the tail of 
the speed distribution. However in these cases the expo- 
sure is large even for the Rayleigh- Watson statistic, due 
to the small values of R(> E t h)- 

We have focused on the number of events and exposure 
required for a '5cr' discovery of light WIMPs with a di- 
rectional detector. Significant experimental support for 
light WIMPs would be obtained even with a lower sig- 
nificance signal. For 95% and 99.713% confidence levels 
(the later corresponding to 3<r) the number of events, and 
hence the exposure, required with the Rayleigh- Watson 



statistic is smaller by a factor of between 0.40 — 0.91 and 
0.69 — 0.95 respectively. The factor by which the number 
of events must be increased to increase the significance of 
the rejection of isotropy is smallest when the experiment 
is only sensitive to the most anisotropic events coming 
from the tail of the speed distribution. However in these 
cases the exposure required even for a low significance 
detection is large. 

The final question is 'How achievable are these ex- 
posures and energy thresholds by current and near- 
future detectors'? A typical current detector consist- 
ing of a 1 m 3 TPC filled to 75 Torr with CF 4 or CS 2 
could, in r oug hly a year, achieve an exposure of order 
10 3 kgday [291 ■ A 10 3 kgday exposure with CF 4 or CS 2 
would (provided that recoils are measured in 3d with 
good, < 10°, angular resolution) be capable of detecting 
m x = lOkeV WIMPs with an energy threshold of 20 keV 
or lower. Lighter WIMPs would require a lower energy 
threshold, potentially as low as 5keV for m x = 5 GeV. 
With a 3 He target a low, ~ 5 keV, energy threshold would 
be required, even for m x — 10 GeV. Measuring the di- 
rections of low energy nuclear recoils is a major exper- 
imental challenge, and these energy threshold are lower 
than those which have been used in the analysis of data 
from current, prototype, detectors [29j]. One of the fo- 
cuses of the R&D for future generation experiments is to 
reduce the energy threshold (see e.g. [59j). The MIMAC 
experiment has, using micromegas readout, detected 5 
keV F recoils [6fj] , while DRIFT-II is sensitive to nuclear 
recoils down to sub 5 keV energies [6l|. Simulations of 
the MIMAC detector indicate that the directional energy 
threshold will lie below 20 keV [H 
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IV. SUMMARY 

The event rate excess and annual modulations ob- 
served by various direct detection experiments may be 
due to light, 5 - 10 GeV WIMPs. We have investi- 
gated whether near future directional detection experi- 
ments will be able to test this possibility, by detecting 
the anisotropy of the nuclear recoils. We find, using the 
Rayleigh- Watson statistic, that an ideal directional de- 
tector (capable of measuring the directions of the recoils, 
including their senses, in 3d with good angular resolu- 
tion) would require between 7 and 58 events to detect 
the anisotropy at 5c The number of events required de- 
pends on the target nuclei mass, energy threshold, WIMP 
mass, and crucially the details of the high speed tail of the 
WIMP speed distribution. It is smallest for cases where 
the experiment is only sensitive to the highly anisotropic 
recoils from high speed WIMPs in the tail of the speed 
distribution. If the detector is not capable of measuring 
the senses of the recoils we find, using the Bingham statis- 
tic, that the number of events required ranges between 9 
and more than 190. The increase in the number of events 
required with axial data is smallest for the cases which 
are only sensitive to high speed WIMPs due to the higher 
degree of anisotropy in the nuclear recoil distributions. 

In terms of the detection potential the key quantity is 
the exposure necessary to detect the required number of 
events, which is inversely proportional to the event rate 
above threshold. For the configurations we have consid- 
ered ( 3 He, C, F and S targets, 5 - 10 GeV WIMPs and 
energy threshold between 5 and 20keV) the event rate 
above threshold varies by more than five orders of mag- 
nitude, and in some cases the minimum speed required 
to cause a recoil above threshold exceeds the maximum 
WIMP speed in the lab and the event rate is zero. For 
the cases which are only sensitive to high speed WIMPs, 
where the number of events required was smallest, the 
event rate above threshold is small and hence very large 
exposures would be required. The shape of the high 
speed tail of the WIMP distribution is not well known 
and this leads to large uncertainties in the event rate 
expected in experiments which are only sensitive to the 
high speed tail, c.f. Refs. [U, H S3 Ell ■ We also em " 
phasize that including, and averaging over, the Earth's 
orbit is essential for an accurate calculation of the event 
rate in these cases. 

We find that a future CF 4 or CS2 detector with an 
energy threshold of 20keV or lower, which can mea- 
sure recoil directions and senses in 3d with good an- 
gular resolution, would be capable of detecting WIMPs 
with (m x ,(T p ) = (10 GeV, 10 -4 pb) with an exposure 
of 10 3 kg day. Detecting lighter WIMPs would require 
a lower energy threshold. With a 3 Hc target a low, 
~ 5keV, energy threshold would be required, even for 
WIMP masses at the upper end of the mass range con- 
sidered. In summary, we conclude that directional detec- 
tion experiments may be able to detect light WIMPs, but 
this depends quite sensitively on both the experimental 



configuration (target nuclei mass and energy threshold) 
and the unknown WIMP mass and velocity distribution. 
The directional energy thresholds required to detect light 
WIMPs are below those used in analyses of data from 
current directional detectors, however it has been demon- 
strated that directional detectors can detect HH and 
measure the directions [34j of low energy nuclear recoils 
and R&D is underway to realise lower energy thresh- 
olds m. 

We also note that directional experiments, with heavy 
targets, could also test inelastic dark matter as the ex- 
planation of the direct detection anomalies [62[ . 
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Appendix A: Rayleigh Watson statistic 

The (modified) Rayleigh- Watson statistic, W* is the 
simplest coordinate independent statistic for detecting 
anisotropy in vectorial data. It is related to the Rayleigh 
statistic, 1Z, which for a sample of N unit vectors 5% is 
given by 



n = 



N 



(Al) 



i.e. the modulus of the sum of vectors. For an isotropic 
data set TZ should be zero, modulo statistical fluctua- 
tions. For anisotropic data the value of TZ becomes larger 
as the degree of anisotropy increases. 

The modified Rayleigh- Watson statistic, W*, defined 

as EMM 



w* = ( 1 - — 1 w + — !— 

2N J ION 



W 2 , 



(A2) 



where W is the (unmodified) Rayleigh- Watson statistic 

3 



W = —TV 

N 



(A3) 



The modified statistic W* has the advantage of approach- 
ing its large N asymptotic distribution for smaller N than 
the unmodified statistic. For isotropically distributed 
vectors, W* is asymptotically distributed as x§ [63l. Ia"ij|. 
The difference between X3 an( A the true distribution of 
W* for isotropic vectors in the large W* tail of the dis- 
tribution is less than 2% for N > 30 (22J. For smaller 
N the xl distribution significantly underestimates the 
true probability distribution and therefore, as in Ref. 22 1 
we calculate the probability distribution from the exact 
probability distribution of TZ, as described in Ref. 
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Appendix B: Bingham statistic 

The Rayleigh- Watson statistic can not be used with 
axial data, as it is not sensitive to distributions which 
are symmetric with respect to the centre of the sphere. 
For axial data the Bingham statistic B* which is based 
on the scatter matrix of the data, T, can be used. This 
matrix is defined as EH, E3] 

2 / x^Xi x^yi X' t z^ \ 
T=— ViX, yiyi y l z i , (Bl) 

i—l \ Z^X^ Zi]Ji ZiZi J 

where (xj, Zj) are the components of the i-th vector or 
axis. This matrix is real and symmetric with unit trace, 



so that that the sum of its eigenvalues (k = 1, 2, 3) is 
unity, and for an isotropic distribution all three eigenval- 
ues should, modulo statistical fluctuations, be equal to 
1/3. The Bingham statistic, B, 

fc=l v 7 

measures the deviation of the eigenvalues from the 
value of 1/3 expected for an isotropic distribution. For 
isotropically distributed vectors/axes B is asymptotically 
distributed as x|- Since T is symmetric under a sign 
change of x, the Bingham statistic can be used for axial 
data as well as vectors. 
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